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. Abstract 

o 

I The similarity theory of stellar structure is used to study the prop- 

■ erties of very massive stars when one can neglect all the sources of 

, opacity except for the Thomson scattering. The dimensionless inter- 

nal structure of such stars is practically independent of the energy 
Q H i generation law. It is shown that the mass-luminosity relation can be 

q approximated by an analytical expression which is virtually universal 

\ regarding the chemical composition and energy generation law. A de- 

tailed comparison with the Eddington standard model is given. The 
application of the obtained results to the observations of massive stars 



is briefly discussed. 

1 Introduction 



The similarity theory of stellar structure was a powerful tool to understand 
the properties of stellar models during the pre-computer era of their study. 
Suffice it to mention a theoretical explanation of the mass-luminosity and 
mass-radius relations (Biermann, 1931; Stromgren, 1936; Sedov, 1959). For 
years the similarity theory remains to be useful for the interpretation of 
various aspects of stellar structure (Schwarzschild, 1958; Chiu, 1968; Cox, 
Guili, 1968; Dibai, Kaplan, 1976; Kippenhahn, Weigert, 1990). 

Here, first we describe the similarity theory of stellar structure as a 
boundary- value problem formulated by Imshennik and Nadyozhin (1968). 
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Then we discuss the structure of chemically homogeneous stars of so large 
a mass that the opacity can be considered as being due to the Thomson 
scattering alone. Such an approximation is believed to be adequate for still 
hypothetical Pop III stars and for a number of observed luminous stars - 
e.g., massive O-stars, Wolf-Rayet stars, and some specific stars like rj Car, 
and the Pistol star. A special consideration is given to the comparison with 
the Eddington standard model. 



2 Similarity theory of stellar models 

Let us consider a spherical star in hydrostatic and thermal equilibrium. If 
the pressure P is contributed by perfect gas P g and black body radiation P T 
and the rate of energy generation e and opacity k are the power functions of 
the temperature T and density p then the structure of the star is described 
by the following set of differential equations: 



(!) 

= 4vrr 2 p , (2) 
4irr 2 pe, (3) 



dP _ Gm 

dr ^ r 
dm 

dr 
dh 

dr 

{pT Gm yj ^ yj yj ( d log T 
P r z \d\ogP J s u\ 

pTGmyj V7 ^vy Y7 3kLP 

— — v r , v r ^v A , v r - 167rcaT 4 Gm , 

P = P g + Pr = —-pT + \aT\ (5) 
m u p 3 

o P z fl-i-14. „ am » T3 v 2(4-3/3) 

/€ = /€ p a T^, £ = £ P 5 T\ (7) 

where /3 is the ratio of the gas pressure to the total one; G, k, and a are the 
gravitational, Boltzmann, and radiation density constants, respectively; m u 
is the atomic mass unit and c is the speed of light. The mean molecular mass 
p and the coefficients k and e depend on composition. 
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The above equations should meet the boundary conditions: 



Center r = : m = , L = , (8) 

Surface r = R : p = 0, p = 0, m = M. (9) 

The conditions (jHJ) tell that there is neither point mass nor point source of 
energy in the center of the star. The conditions (JSJ) at stellar surface state 
that both the pressure and density vanish there and the mass m must be 
equal to the total mass M specified for the star. The stellar radius R is to 
be obtained as a result of the solution of Eqs. (H}-©, i.e. R is an eigenvalue 
of the problem. Simultaneously, the solution gives the luminosity of the star 
L = L(R). 

Assuming that the chemical composition and, consequently, fi, k , and e 
are constant throughout the star, one can reduce the above equations to a 
dimensionless form by measuring the physical quantities in the following set 
of units (Schwarzschild, 1958): 



R, m -> M , L -> L 



o 



(10) 



m u GM „ n GM 2 n M _ 

Introducing the following dimensionless variables 

x = r/R, q = m/M, l = L/L , p = P/P , v = p/p , t = T/T , (11) 
one can rewrite Eqs. (0)-© in the form (Imshennik, Nadyozhin, 1968): 

(12) 
(13) 
(14) 
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where all the parameters are gathered in the three dimensionless constants: 



C 



i 



(4tt) 5 \ k J 
3(47r)- Q f_k 




,' 7 M 1+<5+ ' 7 



L R 3S+r > ' 



(17) 



C 2 



B 



64tt 2 ac \Gm 
Aira ( Gm u \ 4 



3G V k 




(19) 



(18) 



The dimensionless boundary conditions become: 



Center x = : q = , / = 0, 

Surface x = l: p = , t = , / = 1 , g = 1 . 



(20) 
(21) 



Thus, we have now six boundary conditions for four first order differen- 
tial equations ([T2|) - (jl5)) . This means that for every given B (or /i 2 M) the 
constants C\ and C 2 can have only quite certain values (eigenvalues) in order 
that the solution of four differential equations would satisfy to six boundary 
conditions. The solution itself as well as C\ and C 2 depend only on /i 2 M and 
the exponents a, v, 5, r\: Ci j2 = Ci j2 ( / u 2 M, a, v, 5, rf). 

Equations (fT2j) - (pj) can be solved by standard methods. For the calcu- 
lations described in the next section, we have used the following approach. 
Having chosen some trial values of C\ and C 2 and integrating Eqs. ()12|) - 
f!16|) numerically from the surface [x — 1, boundary conditions (J2"T|) ] inward 
to some value x = Xf (0 < Xf < 1), one obtains a set of four quantities 
{p(xf),q(xf),l(xf),t(xf)} as functions of C\ and C 2 . Then integrating these 
equations again from the center [x = 0, boundary conditions (J2(J)) ] with the 
same C\ and C 2 and trial values p(0) = p c and £(0) = t c outward to Xf one 
obtains a similar set of quantities that depends now on C%, C 2 , p c , and t c . 
Since p(x), q(x), l(x), and t(x) must be continuous throughout the star, both 
the sets have to coincide at x = Xf. One can organize an iterative process 
for four unknown parameters Ci, C 2 , p c , and t c to find such their values that 
would ensure the desired continuity of four variables p, q, I, and t at x — Xf. 
Of course, the result does not depend on Xf but the convergence domain of 
the iterations does. 

As soon as C\ and C 2 are found, one can easily derive the luminosity Lq 
and radius R from Eqs. (fTTj) and (|TH|) which represent the mass-radius and 
mass-luminosity relation, respectively. 
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3 Structure of very massive stars 



In this section, the above similarity theory is employed to describe the struc- 
ture of very massive stars. In this case the opacity is dominated by the 
Thomson scattering and one can assume a = 0, v = 0, and k = no = 
0.2(1 + X) cm 2 /g (X is the mass fraction of hydrogen). 

The mass-luminosity relation [Eq. 1)18)1] can be converted into the form 



647T 2 ac ( Gm u 
3 kq 



i 



//M 3 C 2 (/z 2 M,5,7/). (22) 



All the models with physically reasonable values of the energy generation 
temperature exponent (i.e., rj ^ 4) considered below have convective cores 
and radiative outer envelopes. Using Eqs. (J12j) and ()15)) and boundary con- 
ditions (J2U), one can find that C 2 is simply connected with f3 s (ii 2 M,5,r}) - 
the value of the parameter f3 at stellar surface: 

l-(3 s = ABC 2 . (23) 

Now the mass-luminosity relation can be reduced to 

L = L Ed (l - (3 S ) , (24) 

where the Eddington critical luminosity Lgd is given by 

_ _ AttcGM 6.483 x 10 4 M _ , nr . 

LEd = = i+x Mo L& ■ (25) 

Using the Runge-Kutta method with automatic control of accuracy for 
solving Eqs. ()12j l -()16) l from stellar surface down to Xf and from center up to 
Xf, we have calculated a large number of models for a wide interval of the 
parameter /i 2 M (0 < /i 2 M < 4, 000 M ). Typically, the values x f pa 0.1-0.3 
are the best for the iterations to converge. The convergence domain, however, 
turned out to be rather narrow at least for the Newton-Raphson iteration 
scheme used in our calculations. So, when calculating a sequence of such 
models one has to change fi 2 M at most by a few percents in order to ensure 
the convergence. 

Table^presents the most important properties of several selected models. 
Note that for the sake of a better perception we use hereafter for dimension- 
less quantities the same designations as for dimensional ones. The first row 
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Table 1: Structural properties of selected models 
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* CNO-cycle (5 = 1, r) = 16), t 3a-reaction (S = 2, 77 = 30), * pp-chain (5 = 1, 77 = 4) 



of Table Ogives /i 2 M with M measured in Mq . The next three rows are the 
Ci- values for three modes of energy generation: the CNO-cycle, 3a-reaction, 
and pp-chain. The fifth row gives the C^-values. The next six rows relate to 
the central dimensionless values of density p c , pressure P c , temperature T c , 
radiative temperature gradient V rc , gravitational potential <p c (in the units 
GM/R), and ratio f3 c of the gas pressure to the total one. The next two rows 
contain /? conv and (3 a at the convective core boundary and at stellar surface, 
respectively. Finally, the last six rows present the dimensionless values of ra- 
dius, mass, luminosity, density, pressure, and temperature at the convective 
core boundary. 

The calculated sequences of one-parametric (depending on \j?M) dimen- 
sionless models have so large convective cores that the major part of ther- 
monuclear energy turns out to be released within the core. Therefore, L(r) 
becomes almost equal to the luminosity Lq everywhere outside the convec- 
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tive core. For such a strong dependence of the energy generation on tem- 
perature as in the CNO-cycle (rj = 16) and in 3a-reaction (rj = 30), one 
can with a high accuracy assume l(x) = 1 in Eq. (JT5j) . Consequently, Eq. 
(fT3jl for l(x) proves to be disentangled from other equations (JUJ), JUJ), (fTSj) . 
and ()16jl. Now one can calculate the overall stellar structure by solving a 
truncated eigenvalue problem that pays no attention to the law of energy 
generation and considers the constant as a single eigenvalue parameter 
involved. In particular, the resulting ^-eigenvalue, and consequently L 
[Eq. (jUJ)], as well as the dimensionless stellar structure do not depend on 
the exponents 5 and rj. Such a truncated model is actually the point-source 
model with the Thomson opacity and allowance for the radiation pressure 
that was first calculated by Henrich (1943) within a limited range of stellar 
mass (0 ^ fi 2 M ^ 119). Our numerical results are in excellent agreement 
with his ~ 1% accurate calculations. 

When C*2 and the dimensionless functions cr(x) and t(x) are already known 
from the solution of the truncated problem one can find the constant C\ and 
spatial distribution of the dimensionless luminosity l(x): 



Ci 



~ />1 -l 1 px 

/ x 2 a 1+s t r, dx , l{x) =d / x 2 <r 1+s f>dx. (26) 
Jo J Jo 



If the energy generation is an arbitrary (but still strongly varying with tem- 
perature) function: e = e(p,T), then instead of C\ one has to consider the 
stellar radius R as an eigenvalue that can be found from the equation 



M I x 2 a(x) e 



M m u GM 
;cr{x), \i— —t\x) 



AttR 3 v n n k R 



dx , (27) 



where for given M and composition the luminosity L is determined by the 
C 2 -value obtained for the truncated model. 

The structural properties of the models presented in Table [T] correspond 
to three modes of energy generation: the CNO-cycle, 3a-reaction, and pp- 
chain. For the CNO-cycle and 3a-reaction, the dimensionless luminosity 
-^conv at the convective core boundary is equal to 1.0000 for all the values 
of /j, 2 M. The same is true for the pp-chain when fi 2 M > 30. All other 
parameters in Tables [T] and El except for C\ and V rc , prove to be the same 
for the CNO-cycle and 3a-reaction, and at /x 2 M > 30 for the pp-chain as 
well. Practically, one can consider only the interval ^ fi 2 M < 10 as 
exhibiting slight differences for the pp-chain shown in Figures EUSlEl below. 
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r 



Figure 1: Density versus radius in the units defined by Eqs. (jlOJ) for /i 2 M = 
(dashed line), 4000 (full line) and for a polytropic model {n = 3, dotted line). 
The boundaries of the convective core are indicated for fi 2 M = and 4000. 

This concerns also the mass-luminosity relation. More specifically, for t] = 4: 
log 10 C2 = —3.516 and —3.285 at /i 2 M = 10 and 0, respectively. 

The constant C% and central value of the logarithmic radiative tempera- 
ture gradient V rc = CiC 2 p c crf t^~ 4 strongly depend on the energy generation 
law. The V rc values for the CNO-cycle in Table Q demonstrate that even for 
low masses (/x 2 M < 1) when one can neglect the radiation pressure, there is 
still a good excess of V rc over Va (0.25 ^ Va ^ 0.4) to drive the adiabatic 
convection along stellar core. For the 3a-reaction, the excess is much larger: 
V rc = 11-05 for fi 2 M = 0. However for the pp-chain we have V rc = 0.594 - 
the value exceeding the adiabatic gradient only by a factor of 1.5 (Va c = 0.4 
for fi 2 M « 0). This happens because the exponent rj = 4 is not very far 
from a critical value of rj w 2 — 3 necessary to ensure the existence of the 
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Figure 2: Temperature versus radius in the units denned by Eqs. fllOl) for 
fi 2 M = and 4000. The Eddington model is shown by the dotted lines. 

convective core (Cowling, 1934; Naur, Osterbrock, 1953). 

Figures [T] and |2 show the density and temperature versus radius for two 
extreme values of /x 2 M in the case of the CNO-cycle and 3a-reaction. With a 
high accuracy these two modes of energy generation result in only one model 
for every given /z 2 M (see the discussion above). 

It is instructive to compare the results with the Eddington standard model 
(Eddington, 1926) that at constant opacity corresponds to a constant rate 
of energy generation: e = e$ (r) = 0, 5 = 0). In this case G\ — 1 as it 
immediately follows from Eq. (|26p. Then, it is easy to make sure that (3 is 
to be constant inside such a model. As a result, the pressure turns out to be 
connected with density along the radius by a power law: 

P = Kp^ 3 , K= — 

m u /i 
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3k 1-/3 



i 



am n fi j3 A 



const . 



(28) 



Thus, the standard model is nothing else than a polytropic gas sphere of 



. i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i . 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

io gl0 <yM) 

Figure 3: The ratios of the gas pressure to the total one, /3, at the center 
P c , the convective core boundary /3 conv , and the surface (3 a as functions of 
/x 2 M (M is in Mq ). A dotted curve corresponds to the Eddington standard 
model [Eq. flU]- 

index n = 3. The value of j3 proves to be uniquely connected with the mass 
M by the following equation (Chandrasekhar, 1939): 

^ = 0.01607G 3 (^) 4 (^M) 2 = 2.994XKT 3 {f^f ■ (29) 

The constant Ci as a function of /i 2 M and the mass-luminosity relation for 
the standard model are determined by Eqs. (J23j) and (|24*j) with j3 s being 
substituted for (3 from Eq. (J25j) . It is easy to verify that the Eddington 
standard model has no convective core at all — it is convectively stable for 
any fi 2 M. 
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According to Fig. |3J (3 increases from the stellar center up to the surface: 
(3 C < Pconv < As j the Eddington (3 always remaining between f3 c and (3 com . 




Figure 4: The dimensionless central pressure P c , density p c (right vertical 
axis), and temperature T c (left vertical axis) as functions of p 2 M (M is m 
M Q ). 

Figure |U shows the central pressure P c , density p c and temperature T c as 
functions of p 2 M. At p 2 M < 10, these quantities for the pp-chain (dotted 
lines) only slightly differ from those for the CNO-cycle and 3a-reaction (solid 
lines). The broken line represents T c for the Eddington model for which P c 
and p c are independent of p 2 M: P c = 138.9, p c = 162.5 (polytrope n — 3 !). 

Figure |3] shows how the dimensionless radius of the convective core r conv , 
the value of (3 at the convective core boundary /? conv , and the mass above the 
convective core (1 — m conv ) depend on p?M. At p 2 M > 40, the convective 
core contains more than 85% of stellar mass. One can approximate m conv (for 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



log 10 (^i 2 M) 



Figure 5: The convective core boundary characteristics as functions of fi 2 M 
for the CNO-cycle and 3a-reaction (solid lines) and pp-chain (dotted lines). 
The values for the CNO-cycle and 3a-reaction at /i 2 M = are shown explic- 
itly 



all the three energy generation modes!) by the following asymptotic relation 
(a straight broken line in Fig.|3J): 

m conv = l - 9.075 (jj 2 M) ~ L095 , /i 2 M > 40 . (30) 



Since the constant C\ strongly depends on the energy generation mode 
(Table HJ), it is useful to introduce a modified constant C\ m : 

C lm (/i 2 M) = (/i 2 M)^ [(1 - (3 S ) Ci]-^ , (M in Mq ) . (31) 
Now, using Eq. (J24j) for L we can obtain the following expression for stellar 
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Figure 6: The modified constant G\ m as a function of jj?M for the CNO-cycle, 
3a-reaction and pp-chain. The open circles and dashed line show polynomial 
fits. 



radius function of /i 2 M and composition: 

\ogR = \ogC lm + hgD + logS 1 , (R in R Q 

where 

logD 

oo -|- lj 



' l O g[0.2(l + X)e]-§^log/i. 



log S 



35 + rj 
V 



log 



35 + 77 
log(4vr) 



(32) 



35 + rj \ k J 35 + r] 

3o + ?7 VcG/ 35 + r] 
X 2 (pp-chain), XXqno (CNO-cycle), Y 3 (3a-reaction) . 
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Here, the dependences of Eq on the mass fractions of hydrogen X, the CNO- 
isotopes X CNO , and helium Y are shown explicitly (e = £00 0- F° r large 
values of rj, Ci m becomes independent of the energy generation mode - 
compare the CNO (i] = 16) and 3a (r/ = 30) curves in Fig. El Both the Ci m 
can be approximated by a single polynomial shown by open circles in Fig. |3 

y = 0.0720052 + 0.782547a; -0. 120854a; 2 + 0.0295923a; 3 -0.0030574a; 4 , (33) 

where y = log 10 Ci m and x = log 10 (fi 2 M). A fit for the pp-chain (r/ = 4) is 
shown in Fig. El by a dashed line that gives a good accuracy for /i 2 M > 10. 
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Figure 7: The normalized mass-luminosity relation [Eq. ()24|) ]. Dotted curves 
indicate asymptotics for small and large /z 2 M-values. The Eddington model 
is shown by a broken line. 
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Figure 8: The mass-luminosity relation (M in Mq , see text). 



Figure shows luminosity L in terms of L-^d as a function of /i 2 M [solid 
curve specified by Eq. (124 j) ]. The curve virtually holds for all the three energy 
generation modes — a small difference for the pp-chain at /x 2 M < 10 is 
indistinguishable on the scale of the figure. For given /i 2 M, the Eddington 
model is always over-luminous (by a factor of two at fi 2 M < 10). 

For the sake of practical use, it is worth to rewrite the mass-luminosity 
relation in terms of variable A: 



A 



M Ln 



// 



Mq L Ed 



AncGM, 







1.5426 x 10~ 5 /x 2 (l + X) 



(34) 



Connecting the asymptotics for small and large /x 2 M-values by a cubic spline 
that ensures the continuity of the function and its first derivatives (open 
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Table 2: Integral properties of selected models 
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0.0139 


0.0439 


0.0908 


0.1247 


0.1309 


0.1271 


0.1214 




t 8 


2.438 


2.517 


2.580 


2.624 


2.649 


2.665 


2.677 


2.723 


tconv 


0.269 


0.433 


0.619 


0.843 


1.044 


1.251 


1.472 




r c 


1.441 


1.426 


1.475 


1.664 


1.896 


2.132 


2.329 


2.587 


<7> 


5/3 


1.532 


1.453 


1.401 


1.3733 


1.3555 


1.3445 


1.345 


UJ 


2.804 


2.109 


1.629 


1.288 


1.058 


0.840 


0.627 


0.667 



* E g attains a maximum —1.2005 at /i 2 M = 18.5 



/ attains a maximum 0.1636 at fj, 2 M = 21 



circles in Fig. |HJ) one obtains the following analytical approximation: 

lg(/i 2 M) = 0.9347 + |lgA, (IgA < -1.7), 

lg(/x 2 M) = 1.015209 + 0.449843 lg A + 0.0534969 lg 2 A (35) 
+ 0.00754094 lg 3 A, (-1.7 < lg A < 1.7), 

/i 2 M = 10.125 (l + 0.0987654A + Vl + 0.1975A) , (IgA ^ 1.7). 

Equation ()35j) allows to estimate M when the luminosity Lq and the repre- 
sentatives of composition (/z and X) are known. On the contrary, in order to 
estimate Lq for given M and composition one can either solve Eq. (|33j) for A 
or use the following accurate approximation: 

A = 0.00157 (zi 2 M) 3 , (/i 2 M *C 2.4) , 

lg A = -2.907029 + 3.552793 lg (/x 2 M) - 0.7717945 lg 2 (/i 2 M) (36) 
+ 0.078623 lg 3 (/i 2 M) , (2.4 < /i 2 M < 100) , 

A = ji 2 M (1 - -4=\ , (H 2 M > 100) . 

Concluding this section we present in Table |2] a number of integral prop- 
erties of stellar models such as gravitational E g and thermal Et energies (in 
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log 10 (|a 2 M) 

Figure 9: Dimensionless gravitational E g and thermal Et energies, central 
moment of inertia /, and gravitational binding energy E\> versus fi 2 M. 

units GM 2 /R), central moment of inertia of the whole star / = J Q M r 2 dm 
and that of the convective core 7 conv (in units MR 2 ), the time sound takes 
to propagate from the center up to the surface t s = J Q fl — ^== and to the 

convective core boundary t conv (in units \J R 3 / (GM) ), column density of the 

star t c = jf^pdr (in units M/R 2 ), average adiabatic index (7) = ^fw^ttt~ ^ 

and an estimate of fundamental angular frequency cu = a/ (3(7) — 4)|_E g |/J 
of radial pulsations (in units a/ GM/R 3 ). 

The last column in Table El presents the properties of the Eddington 
standard model at fi 2 M = 4000. Note that for the standard model the 
dimensionless E g , I and r c are determined by the polytrope n = 3 structure 
and do not depend on /z 2 M. 
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Figure 121 shows E g , E?, I, and gravitational binding energy E h = —E tot = 
— (E g + E T ) as functions of fi 2 M. 



4 Comparison with detailed models 

In order to demonstrate potentialities of the similarity theory we compare 
our results with detailed models of massive main sequence stars calculated 
by Schaller et al. (1992) and the models of helium and carbon-oxygen stars 
(Wolf-Rayet stars) studied by Langer (1989) and Deinzer and Salpeter (1964). 
Figures ITUHT31 display such a comparison. Solid curves in Fig. [TUlare obtained 
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Figure 10: The mass-luminosity relation (solid lines) in comparison with the 
detailed models of Schaller et al. (1992) (filled and open circles). M and L 
are in solar units. 



with the aid of Eq. (|36|) for two specified by Schaller et al. (1992) compo- 
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sitions: X=0.680, Y=0.300, Z=0.020 (upper curve) and X=0.756, Y=0.243, 
Z=0.001 (lower curve). One can observe an excellent agreement for low 
metallicity (Z = 0.001) and a fairly good coincidence for Z = 0.02. In the 
latter case our models are slightly over-luminous (by 25% at M — 15 Mq). 
This natural result can be interpreted as on average a 25% contribution to 
the opacity coming from sources other than the electron scattering. 
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Figure 11: The mass-radius relation (solid lines) in comparison with the 
detailed models of Schaller et al. (1992) (filled and open circles). M and R 
are in solar units. 

The mass-radius relation for the same models is shown in Fig. ^2 The de- 
tailed models have systematically larger radii (typically by ~ 5% for Z=0.02) 
as compared with our models (solid lines). This can be explained by a no- 
ticeable increase in opacity in stellar envelope owing to the absorption in 
atomic spectral lines (Imshennik, Nadyozhin, 1967). The solid lines were 
calculated with the aid of Eqs. and (J3l3|) assuming the CNO-cycle as the 
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main energy source with the energy generation law taken from Caughlan, 
Fowler (1988) and approximated by 

£ C no = 8.43xl0- 3 X CNO XpT 7 16 ergg'V 1 , (2.2 < T 7 ^ 3.6) , (37) 

where T 7 = T/10 7 K and p is in gem -3 . Within the indicated limits for T 7 , 
the accuracy of Eq. (}3*Tj) is better than 10%. 

1.0 

0.9 

1 0.8 
o 

E 

0.7 
0.6 
0.5 
0.4 

1.1 1.3 1.5 1.7 1.9 2.1 

logio M 

Figure 12: Dimensionless mass of the convective core m conv and central value 
of P (solid and dashed lines) in comparison with Schaller et al. (1992) (filled 
and open circles and squares). Solid and dashed lines are for Z = 0.001 and 
Z = 0.02, respectively; M is in Mq . 

Similar approximation for the 3a-reaction is given by the following power 

law 

e 3a = 4.95 x 10~ 38 Y 3 p 2 T 30 erg g'V 1 , (T 7 > 8) . (38) 
The accuracy of the approximation is better than 10% at T 7 > 8. 
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The structural properties of our models are also in a good agreement with 
the detailed models as it is exemplified in Fig. [T2]for central value of (3 and 
the convective core mass fraction m conv . 
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Figure 13: The mass-luminosity relation for the helium and carbon-oxygen 
Wolf-Rayet stars (solid lines) in comparison with the detailed models of 
Langer (1989) (filled and open circles) and Deinzer and Salpeter (1964) (as- 
terisks). M and L are in solar units. 

Figure ITBI shows the mass-luminosity relation for pure helium stars (Y=l) 
calculated by Deinzer and Salpeter (1964) and Langer (1989), and for carbon- 
oxygen Wolf-Rayet stars, with the mass fractions of carbon, oxygen, and he- 
lium Xc=0. 113, Xo=0.867, Y=0.02, respectively (Langer, 1989), in compari- 
son with the mass-luminosity relation for our models determined by Eq. (|35|) 
(solid lines). One has to specify the hydrogen mass fraction X=0 in Eq. ()34|) 
for A and make use of the following expression for the mean molecular mass 
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fi: 

^ 3 6 Y + 28X C + 27X - (Y + X c + X = l). (39) 

Our results are in a very good agreement with the detailed models even 
for such a small mass as 1 Mq. 



5 Discussion and conclusions 

It is interesting to apply the mass-luminosity relation given by Eqs. (I35j) 
and ()36j) to the very massive stars observed in the Milky Way, in Magellanic 
Clouds and in a number of nearby resolved galaxies. There exist compre- 
hensive observational data for several dozens of such stars (see, for instance, 
Figer et al., 1998; Puis et a!., 1996; Humphreys, Davidson, 1994; Sandage, 
Tammann, 1974). A detailed analysis of the efficiency of the similarity the- 
ory in deriving the properties of such stars from observations is in need of 
a special paper. Here as an example we consider the Pistol star studied in 
detail by Figer et al. (1998). Assuming that the star was initially of so- 
lar composition (X=0.707, Y=0.274, Z=0.019; Anders, Grevesse, 1989) and 
had the luminosity L = lO 6 7±a5 L (Naj arro, Figer, 1999), one can esti- 
mate from Eq. (|35j) its initial mass to be M = 116, 246, and 595 Mq for 
L = 10 6 ' 2 , 10 6 ' 7 , and 10 7 ' 2 Lq, respectively. The corresponding initial radii 
derived from Eqs. (J3*2j) and turn out to be R = 14, 22, and 36 Rq for 
the CNO energy generation rate given by Eq. IJ37J) . The range of masses 
(116 — 595) Mq corresponds to the range (44 — 230) in the parameter /i 2 M 
since fi = 0.618 for the solar composition. Such stars have luminosity of 
(0.4 — 0.7)LEd (Fig. Ej) and very large convective cores: m conv = 0.86 — 0.976 
[Eq. ()30|) or Fig. Ej. Estimating the Pistol star initial mass we have implied 
that its luminosity has not changed appreciably during the evolution. The 
star seems to be in a state close to the exhaustion of hydrogen in the con- 
vective core — i.e., it is about to leave the Main Sequence. According to 
Schaller et al. (1992), an increment in the luminosity ALo within the Main 
Sequence strip decreases with increasing mass M. Since for M = 120 Mq 
(solar metallicity) ALo ~ 25%, initially the Pistol star might have been by 
~ 10% less massive than the estimates stated above. 

The properties of stellar structure obtained here are in a fair agreement 
also with the detailed models of the very massive, (100 — 250) Mq, initially 
zero-metallicity Pop III stars calculated in the UCSC astrophysical group 
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(Woosley, 2005), especially when the stars are definitely settled at their Main 
Sequence in the state of thermal equilibrium. 

The similarity theory of the structure of very massive stars formulated 
here allows to obtain the following important results: 

1. A simple approximation formula for the mass-luminosity relation, given 
by Eqs. (JM}, (13*5]) . and (JHHJ), is valid for different chemical compositions 
and virtually does not depend on the energy generation mode (pp- 
chain, CNO-cycle, or 3a-reaction). 

2. The overall structure of very massive stars, described in dimensionless 
variables listed in Eq. ([11)1 . depends only on the parameter /j?M where 
fi is the mean molecular mass. Practically, this structure turns out to 
be independent of the energy generation mode — be it the CNO-cycle, 
3a-reaction, or pp-chain. 

3. Although with increasing fj, 2 M the stellar structure approaches that of 
the Eddington standard model, the convergence proves to be rather 
slow. Even for fi 2 M = 4000 there are still noticeable discrepancies 
in some stellar parameters. For instance, the central dimensionless 
pressure and density are respectively by about 20% and 15% lower 
than for the Eddington model (Fig. 0J). 
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